Three-neutrino model analysis of the world's oscillation data 
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A model of neutrino oscillation experiments is constructed. The experiments incorporated are: 
solar neutrinos (Chlorine, Gallium, Super-K, and SNO), reactor neutrinos (Bugey and CHOOZ), 
beam stop neutrinos (LSND decay at rest and decay in flight), and atmospheric neutrinos. Utilizing 
this model and the standard three-neutrino mixing extension of the standard model, the data are 
analyzed. Solutions for the mixing angles and mass-squared differences are found to occur in pairs 
corresponding to the interchange Am 2 2 «-* Am^. Two pairs of solutions are found that reasonably 
reproduce the data, including the LSND data. These solutions are #12 ~ 0.5, #13 « 0.1, 622, ~ 0.7, 
Am.12 ~ 5 x 10 -5 eV 2 and Am 2 3 ~ 0.2 eV 2 or 2.4 eV 2 . Other statistically significant solutions are 
also found which produce negligible oscillations for the LSND experiments. 

PACS numbers: 14.60.-z,14.60.Pq 
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Evidence for neutrino oscillations arises from solar neu- 
trino experiments El 12, E| El B, El El B , neutrinos emit- 
ted from reactors |3, llll Il2l. neutrinos from the 
beam stop of an accelerator |l4| , and neutrinos orig- 
inating from cosmic rays impinging on the atmosphere 
U3 Dii D3 D3 We here model what we believe to be the 
essential physics of each of these experiments. We ana- 
lyze this model utilizing the standard three neutrino mix- 
ing extension of the standard model. We find a number 
of sets of parameters, mixing angles and mass-squared 
differences, which reproduce the data, some of which re- 



produce the entire data set including the LSND experi- 
ments. Previous examinations of three neutrino mixing 
have either excluded the LSND experiments 0,H(J, lim- 
ited the mass-s quared differences & or u™pprox- 
imations and constraints in order to work analytically 

mm mm. 

Neutrino oscillations arise because neutrinos are cre- 
ated in flavor eigenstates, and the flavor eigenstates are 
not equal to the mass eigenstates. The flavor eigenstates 
(labeled by a) are related to the mass eigenstates (labeled 
by k) through a unitary matrix U a k, 



O 
=5 



X 



C12C13 

U a k — > I -S12C23 - C12S23S13 
S12C23 - C12C23S13 



where c a k = cos 8 a k, s a k — sin 6 a k, and 8 a k is real. We 
assume no CP violation. We order the mass eigenstates 
by decreasing mass, and the flavor eigenstates are ordered 
electron, mu, tau. The probability that a neutrino of 
flavor a will be detected a distance L away as a neutrino 
of flavor j3 is given by 



P a ^(L/E) 
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l.27Arrij k L/E, where the units of L/E are 



2 2 



in units of eV 2 . 



For the neutrino flux emitted by the sun we use the 
standard solar model |^. A number of detectors have 
measured the flux of solar neutrinos arriving here on 
Earth. Each has a different acceptance and thus mea- 
sures different energy neutrinos. Each measures a deficit 
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as compared to the flux predicted by the standard so- 
lar model. In order to reproduce the energy dependence 
of the survival rate of electron neutrinos arriving at the 
Earth as seen in the experiments, we must invoke the 
MSW effect Ei The MSW effect arises because 
the neutrinos created in the sun propagate through a 
medium with a significant electron density. The for- 
ward coherent elastic neutrino-electron scattering pro- 
duces an effective change, relative to the mu and tau 
neutrino, in the mass of the electron neutrino given by 
A(r) — \/2G E p(r)/m n , with p(r) the electron density 
at a radius r, G the weak coupling constant, and m n the 
nucleon mass. In the flavor basis, the Hamiltonian then 
becomes 



H mat = UMW+A, 



(3) 



with M. the (diagonal) mass-squared matrix in the mass 
eigenstate basis and A the 3x3 matrix with the interac- 
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tion A(r) as the electron-electron matrix element and ze- 
roes elsewhere. By diagonalizing this Hamiltonian with a 
unitary transformation D ak (r, E) , we define local masses 
and eigenstates as a function or r and E. Care must be 
taken so that D ak (r, E) becomes U a k in the limit of zero 
electron density. In the adiabatic limit, which we use, 
the electron survival probability is 

P^E)=j2D(r,E)l k U 2 ek . (4) 
fc=i 

Neutrinos are produced throughout the sun by vari- 
ous reactions, each with its own energy spectrum. The 
surviving neutrinos are then detected by detectors which 
have a different acceptance for each energy of the neu- 
trino. We model this by taking the survival probability 
for an electron neutrino in a particular experiment to be 
given by 

p ee = £pf / /iW dr / 9j(E) P»"(r, E) dE . 

(5) 

Here, j labels a particular nuclear reaction; we include 
three reactions - pp, 7 Be, and 8 B. The quantity pj x is the 
probability that in a particular experiment the neutrino 
arose from nuclear reaction j. We take these from the 
analysis of Ref. |29j for the solar experiments: chlorine 
3, gallium (GALLEX/ GNO,Sage) % 0,0,0, Super-K 
5j, and SNO ||. The function fj(r) is the probability 
that a neutrino is created by reaction j at a radius r |26j 
of the sun and is integrated from the center of the sun 
to the solar radius Rq. The function gj(E) is the en- 
ergy distribution of the neutrinos emitted in reaction j. 
For 7 Be this is a delta function at 0.88 MeV. The lower 
emission line does not contribute significantly. For the 
pp neutrinos, the energy distribution times the detector 
acceptance is a relatively narrow function of energy; we 
set E to its average. For 8 B neutrinos, we use the en- 
ergy distribution from the standard solar model |26l | and 
numerically perform the integration. 

We check our treatment of solar neutrinos by perform- 
ing a two neutrino analysis of only the solar neutrino 
experiments. We find a minimum for tan 2 9 = 0.43 
with bounds 0.28 < tan 2 9 < 0.68 and a minimum for 
Am 2 = 4.0 x 10~ 5 eV 2 with bounds 1.0 x 10" 5 eV 2 < 
Am 2 < 8.6 x 10~ 5 eV 2 . The analysis from Ref. 
gives tan 2 9 = 0.44 with bounds 0.36 < tan 2 9 < 0.58 
and a minimum for Am 2 = 7.0 x 10 -5 eV 2 with bounds 
5.9 x 10~ 5 eV 2 < Am 2 < 9.7 x 10~ 5 eV 2 . This is sat- 
isfactory for our goal of locating possible solutions that 
are semi-quantitatively correct. Our errors are neces- 
sarily larger than in a thorough and model-independent 
analysis as we do not include all of the data, such as 
the measured neutrino energy spectra. The experimen- 
tal data for the solar experiments which we fit are given 



Experiment 


Measured 


L/E (m/MeV) ) 


Data 


LSND-DAR 




.73 


3.1 ± 1.3 x 10" 3 


LSND-DIF 




.40 


2.6 ± 1.1 x 10~ a 


CHOOZ 


^ee 


300. 


> 0.96 


Bugey 


^ee 


10.3 


> 0.95 


KamLAND 


^ee 


4.1 x 10 4 


.611 ± .094 


Super-K 


T^ee 


2.2 x 10 1U 


.465 ± .094 


SNO 


^ee 


2.2 x 10 lu 


.348 ± .073 


Chlorine 


^ee 


4.0 x 10 1U 


.337 ± .065 


Gallium 


/ ee 


OO. X 1U 


.000 ± .U4o 


Atmospheric 




1.9 x 10 4 


1.13 ±.09 


Atmospheric 


TZ e 


2.6 x 10 3 
o.U 


.85 ± .16 


Atmospheric 




1.6 x 10 4 
18.8 


.73 ± .06 


Atmospheric 




4.3 x 10 4 
5.0 


.62 ± .09 



TABLE I: The experiments, quantity measured, the average 
value of L/E, and experimental data for those quantities in- 
cluded in the model are presented. For the atmospheric data, 
the quantity 1Z is define in the text, and the upper (lower) 
value of L/E is for upward (downward) going neutrinos. Note 
that Super-K plus SNO combined provide a measurement of 

Vee and Vee +Ve„ + Ver ■ 

in Table H| We take values from [2^| which differ slightly 
from the original analysis. 

The reactor ex per iments that we include are Bugey 
0, CHOOZ [Hdll and KamLAND [l^. As there is 
an energy distribution for the neutrinos emitted from a 
reactor, the electron survival probability given by Eq. [21 
must be integrated over this spectrum. For small values 
of L/E, the coherent limit, Eq. [21 remains correct. For 
sufficiently large values of L/E, the incoherent limit, the 
sin 2 (f> osc term averages to 1/2. The transition between 
these regions depends on the details of the energy dis- 
tribution of the source neutrinos. We simplify this by 
using an average value for L/E and by using sin 2 4> osc for 
(j) osc < tt/4 and setting sin 2 cj> osc =1/2 for cj> osc > ir/4. 
KamLAND is unique among these as it measures V ee 
where the others set limits. 

The LSND experiments use neutrinos produced from 
muons created in the LAMPF beam stop. There are two 
experiments. The decay at rest experiment [l3| mea- 
sures the oscillation of an muon antincutrino into an 
electron antincutrino, while the decay in flight experi- 
ment fl4| measures the oscillation of a muon neutrino 
into an electron neutrino. These experiments are unique 
in that they measure the appearance of a different flavor 
neutrino rather than the disappearance of electron neu- 
trinos. We treat the sin 2 <fr osc in Eq. [21 just as we did for 
the reactor neutrinos. 

The Super-Kamiokande experiment 0, [iH [Til Hsf has 
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also measured neutrinos that originate from cosmic rays 
hitting the upper atmosphere. The detector distinguishes 
between e-like (electron and anti-electron) neutrinos and 
/i-like (muon and anti-muon) neutrinos. The rate of e-like 
neutrinos of energy E arriving at the detector from a 
source a distance L away is 



Re(L, E) — Vee{L, E) + n(E)V ef _i{L, E) , 
and for /z-like neutrinos 

R tl (L,E)=V^(L,E) 



(0) 



1 

n(E) 



V efl (L,E), (7) 



where n(E) is the ratio of /i-like neutrinos to e-like neu- 
trinos at the source. We separate the data into neutrinos 
with an energy less than or greater than 1 GeV. An aver- 
age value for the energy is calculated from results given in 

Ref. niiia 

which uses the model of the neutrino fluxes 
from Ref. |3l|. We approximate the energy distribution 
of these neutrinos by fitting a simple "teepee" shaped 
function to the distributions given in Ref. [lfil Il7l | . This 
allows us to do the energy integral analytically. The en- 
ergy averaged values of R a we call R a (r, (E)) The high- 
energy /x-like events are classified as "fully contained" or 
"partially contained" events and each of these arises from 
a different energy spectrum. We combined these as 0.24 
fully contained and 0.76 partially contained [Til ]. 

To remove the dependence on the overall normalization 
of neutrino flux, the ratio of measured fluxes for upward 
going (coming from the far side of the earth) neutrinos 
to downward going (coming from overhead) neutrinos is 
taken. This ratio is 



Tin 



Ra(r up , (E)) 
{rdown, (E)) 



(8) 



We utilize a definition of upward/downward going neu- 
trinos as those with the scattered lepton direction in the 
detector of no more than n/5 radians off-axis. The down- 
ward going neutrinos were assumed to travel 15 km from 
the top of the atmosphere, whereas the upward ones 
travel one earth diameter, 13,000 km. The experiment 
measures the neutrino fluxes as a function of the az- 
imuthal angle. We utilize only the endpoints of these 
measurements. Assuming that if we fit the endpoints, 
the curve between would equally well be fit, we divide 
the error associated with the endpoint by \/5/2 (5/2 = 
the number of experimental points over the number used) 
to more properly weight this experiment with respect to 
the others. The data and the parameters used for all the 
experiments are given in Table [I] 

We fit the mixing angles and the mass-squared differ- 
ences to the data by minimizing chi-squared per degree 
of freedom Xdof- ^ n Table ITT1 we present the parameters 
for the ten best fits. We have also found local minima 
for Xdof near 2.3 and 2.7. Notice that the solutions come 
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.001 
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tt/2 


.63 


1.2 


1.2 x 10" J 


1.4 x 10~ b 
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1.4 x 10~ b 


.11 


10 


2.0 
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TABLE II: The value of Xdof an d the fit parameters for the 
ten best fits. 



and 10) corresponding to the interchange Amf 2 <-* Am^ 
and an appropriate redefinition of the mixing angles. We 
derive this symmetry elsewhere l32j . 

The experimental data and the theoretical results for 
these fits are presented in Table Hnl Note that solutions 1 
through 4 produce non-negligible and reasonable results 
for the LSND experiments while fitting the remainder of 
the experiments. We can see how this comes about by 
first considering the solar experiments. We find a mass- 
squared difference on the order of 10 -5 eV 2 , which is the 
magnitude needed for the MSW effect to produce the en- 
ergy dependence of the solar neutrinos. Quantitatively, 
this is different from the 2x2 case as the MSW effect is 
producing an energy dependence for a case where V ee {oo) 
is about 0.6. In fact, we find that KamLAND is measur- 
ing P ee (oo). 

Considering LSND and the reactor experiments Bugey 
and CHOOZ [llj|, we note that the term in Eq.Hwith 
Am 2 w 10~ 5 eV 2 does not contribute. The coefficients 
of sin 2 <f) osc for V eil for the large mass-squared difference 
terms are each of the order 10 -2 . They are of oppo- 
site sign and equal to each other to about ten percent. 
This produces the LSND results for of order 10~ 3 . 
For V ee these coefficients are of the same sign but are 
individually of the order of 10~ 3 and thus do not give re- 
sults which contradict Bugey or CHOOZ. The coefficient 
of the sin 2 <jf sc term corresponding to the mass-squared 
difference of 10 -5 eV 2 term is about 0.25. This term then 
contributes significantly to the asymptotic value of V ee 
which is important for fitting KamLAND and the solar 
experiments. For solutions 1 and 4, <f> osc for the large 
mass-squared difference terms for LSND experiments are 
in the coherent region, while for solutions 2 and 3 they 
are near tt/2. 

Finally we look at the atmospheric data. We find that 
for solutions 2 and 3, the solutions with the largest mass- 
squared difference, all the values of (jf sc are greater than 
7r/2. This does not give a perfect fit to the atmospheric 
data, but it is sufficiently close that when combined with 



in pairs, (solutions 1 and 4, 2 and 3, 5 and 6, 7 and 8, 9 an excellent fit to the LSND data, a good Xdof results. 
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3.1 ± 1.3 


2.4 


2.8 


2.7 


1.5 


.23 


.13 


.00 


.00 


.04 


.04 


LSND-DIF (xlCT^) 


2.6 ± 1.1 


.74 


2.8 


2.7 


.45 


.07 


.04 


.00 


.00 


.01 


.00 


CHOOZ 


> .96 


.95 


.99 


.99 


.96 


.98 


.99 


1.0 


1.0 


.99 


1.0 


Bugey 


> .95 


.95 


.99 


.99 


.96 


.98 


.99 


1.0 


1.0 


.99 


1.0 


KamLAND 


.611 ± .094 


.58 


.57 


.57 


.55 


.57 


.57 


.60 


.60 


.59 


.59 
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Afi.^ -A- HQ/1 
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.04 


.00 
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.00 
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.00 


.00 
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.348 ± .075 


.33 


.34 


.35 


.34 


.33 


.34 


.35 


.35 


.35 


.35 


Chlorine 


.337 ± .065 


.39 


.39 


.40 


.39 


.38 


.38 


.37 


.37 


.37 


.37 


Gallium 


.550 ± .048 


.54 


.54 


.54 


.52 


.53 


.53 


.49 


.49 


.49 


.49 




1.13 ±.09 


1.17 


1.27 


1.26 


1.15 


1.12 


1.12 


1.12 


1.12 


1.06 


1.06 


n e 


.85 ±.16 


1.04 


1.06 


1.07 


1.05 


1.05 


1.05 


1.01 


1.01 


1.01 


1.01 




.73 ± .06 


.79 


.85 


.86 


.82 


.73 


.73 


.73 


.73 


.73 


.72 




.62 ± .09 


.64 


.69 


.70 


.66 


.62 


.63 


.62 


.62 


.65 


.66 



TABLE III: The experimental results and the predictions of the model for the ten fits given in Table ITT1 



For solutions 1 and 4, the large mass-squared difference 
is smaller in order to better fit the atmospheric data, but 
at a cost to the fit of LSND. 

We also note that our solution 7 corresponds to that 
found in Ref. [2(J, in which LSND was not included. This 
helps give us confidence that our simplified model of the 
experiments is capable of locating possible solutions. 

We have found a set of mixing angles and mass-squared 
differences that, within a model, can produce results that 
fit the world's data, including the LSND experiments. 
The future requires a thorough and model-independent 
analysis to see whether this is actually so. We have also 
found a new symmetry that arises from the interchange 
of the mass-squared differences. 

The authors are grateful for very helpful conversations 
with D. V. Ahluwalia and I. Stancu. This work is sup- 
ported by the U.S. Department of Energy under grant 
No. DE-FG02-96ER40963. 
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